Effective medium theory of left-handed materials 
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We analyze the transmission and reflection data obtained through transfer matrix calculations on 
metamaterials of finite lengths, to determine their effective permittivity e and permeability (i. Our 
study concerns metamaterial structures composed of periodic arrangements of wires, cut-wires, 
split ring resonators (SRRs), closed-SRRs, and both wires and SRRs. We find that the SRRs 
have a strong electric response, equivalent to that of cut-wires, which dominates the behavior of 
left-handed materials (LHM). Analytical expressions for the effective parameters of the different 
structures are given, which can be used to explain the transmission characteristics of LHMs. Of 
particular relevance is the criterion introduced by our studies to identify if an experimental trans- 
mission peak is left- or right-handed. 
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Recently, there have been many studies about meta- 
materials that have a negative refractive index n. These 
materials, called left-handed materials (LHMs), theoreti- 
cally discussed by Veselago [1] , have simultaneously neg- 
ative electrical permittivity e and magnetic permeabil- 
ity fi. Such materials consisting of split ring resonators 
(SRRs) and continuous wires were first introduced by 
Pendry [2,3], who suggested that they can also act as 
perfect lenses [4]. 

Since the original microwave experiment by Smith et 
ai. [5], which first materialized Pendry's proposal, vari- 
ous new samples were prepared [6,7] (composed of SRRs 
and wires) all of which have been shown to exhibit a pass 
band in which it was assumed that e and /i are both neg- 
ative. This assumption was based on measuring indepen- 
dently the transmission, T, of the wires alone, and then 
the T of the SRRs alone. If the peak in the combined 
metamaterial composed of SRRs+wires were in the stop 
bands for the wires alone (which corresponds to negative 
e) and for the SRRs alone (which is thought to corre- 
spond to negative fi) the peak was considered to be left- 
handed (LH). Further support to this interpretation was 
provided by the demonstration that some of these mate- 
rials exhibit negative refraction of electromagnetic waves 
[8] . Subsequent experiments [9] have reaffirmed the prop- 
erty of negative refraction, giving strong support to the 
interpretation that these metamaterials can be correctly 
described by negative permittivity and negative perme- 
ability. However, as we shall show in the present study, 
this is not always the case. The combined system of wires 
and SRRs exhibits synergy of the two components as a 
result of which its effective plasma frequency, ui' is much 



lower than the plasma frequency of the wires, u> p . 

There is also a significant amount of numerical work 
[10-13] in which the complex transmission and reflection 
amplitudes are calculated for a finite length of metama- 
terial. Using these data a retrieval procedure can then be 
applied to obtain the effective permittivity e and perme- 
ability fi, under the assumption that the metamaterial 
can be treated as homogeneous. This procedure con- 
firmed [14,15] that a medium composed of SRRs and 
wires could indeed be characterized by effective e and 
fi whose real parts were both negative over a finite fre- 
quency band, as was the real part of the refractive index 
n. 

In the present paper, we study the transmission of pe- 
riodic systems made up of wires alone, SRRs alone and of 
combined systems of wires and SRRs (LHMs). Through 
a very detailed retrieval scheme [14], the effective e and 
fi of those systems are obtained. It is shown that the 
SRRs have also an electric response, in addition to their 
magnetic response which was first described by Pendry 
[3] . The electric response of the SRRs is demonstrated 
by closing their air gaps, and therefore destroying their 
magnetic response. In fact, it is shown that the electric 
response of the SRRs is identical to that of cut- wires. 
Analytical expressions for the effective e of wires and 
SRRs as well as for the effective \x of SRRs are given. 
Using these analytical expressions one is able to repro- 
duce the low frequency transmission, T, and reflection, 
R, characteristics of LHMs. Even the minor details in T 
and R observed in the simulations can be analytically ex- 
plained. The main power of the present analysis though 
is that it gives an easy criterion to identify if an exper- 



imental transmission peak is LH or right-handed (RH): 
If the closing of the gaps of the SRRs in a given LHM 
structure removes the peak close to the position of the 
SRR dip from the T spectrum, this is implies that the T 
peak is indeed left-handed. This criterion is very valu- 
able in experimental studies, where one can not easily 
obtain the effective e and [x and holds provided the mag- 
netic resonance frequency uj m is well separated from ui' p . 
Our criterion is used experimentally and is found that 
some T peaks that were thought to be LH, turn out to 
be right-handed [16]. 

We use the transfer matrix method to simulate nu- 
merically the transmission properties of the metamateri- 
als. A representative result is shown in Fig. 1, present- 
ing transmission spectra for a usual SRRs+wires meta- 
material and for its isolated constituents. As expected, 
the SRRs system (Fig. 1(a)) shows a well defined T dip 
near ua/c — 0.04 (a is the discretization length, related 
to the unit cell as explained in the caption of Fig. 1, 
and c is the vacuum speed of light), due to its mag- 
netic resonance, which gives a negative effective /i in the 
dip regime, while its effective e is positive in this fre- 
quency range. This interpretation has been verified with 
the retrieval procedure. The lattice of continuous wires 
alone shows no transmission up to the plasma frequency 
ujpd/c — 0.09 (see the dotted line in Fig. 1(b)) due to a 
negative effective e. Combining the SRRs and the wires 
systems we find what is shown in Fig. 1(b) - solid line: 
no transmission at low frequencies (due to negative e), 
a LH transmission peak (where both fi and e are nega- 
tive) that roughly corresponds to the transmission dip of 
the isolated SRRs, followed by a second transmission gap 
(where \i > 0, e < 0) and terminated by a right-handed 
transmission shoulder (at about 0.06). The position of 
the last right-handed shoulder, which does not coincide 
with the oj p of the isolated wires, shows that the SRRs 
contribute also to the electric response of the combined 
system, something that has never been considered pre- 
viously. This SRR contribution leads to a new plasma 
frequency u>' smaller than u> p . The present interpreta- 
tion of the data is supported again by results obtained 
from the retrieval procedure. Thus, combining the res- 
onant behavior in /i(w) of the SRRs with the negative 
behavior of e(u>) due to the wires only is clearly inade- 
quate; the resonant contribution to e(u>) due to the SRRs 
has to be added. 

Furthermore, the interaction of the electric response of 
wires and SRRs has to be taken also into account. Some 
indication that there exists such a non-negligible elec- 
tric interaction between the SRRs and the wires comes 
from the comparison of the T for the in-plane system (i.e. 
wires next to the SRRs, as in Fig. 1(b)), with the T for 
an off-plane system (wires behind SRRs, see Fig. 1(c)- 
solid line). It is found that for the in-plane case the first, 
left-handed transmission peak is much narrower and that 
the oj' p is higher than for the off-plane case. The reduc- 
tion of u>o (the electric resonance frequency of the SRRs, 
see Eq. 3 below) and consequently of the ui' (Fig. 1(c)) 



suggests that there is a significant coupling between the 
electric response of the wires and that of the SRRs in the 
off-plane case. 
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FIG. 1. (a): Frequency dependence of the transmission co- 
efficient of a lattice of SRRs (solid line) and of closed-SRRs 
(dashed line); (b): As in (a) for a metamaterial of SRRs plus 
wires (solid line) and of closed-SRRs plus wires (dashed line), 
with the wires next to the SRRs. The dotted line shows the 
transmission for wires only, (c) As in panel (b) but with the 
wires in the opposite face of the SRRs. The insets show the 
geometry of the unit cell. The relative permittivity for the 
metals is taken to be e m = (— 3 + i5.88)10 5 and for the dielec- 
tric board et, = 12.3. All relative permeabilities are one. The 
unit cell size is 6a x 14a x 14a. a is the discretization length 
and c the light velocity in air. 

Since all the metals, the vacuum and the dielectric 
boards are non-magnetic materials, the only source of 
major magnetic response of the LHM are the circulat- 
ing induced currents in the metallic rings of the SRRs, 
which are forced to oscillate due to the gaps acting like 
capacitors. Therefore, a simple way to study the com- 
bined electric response of the LHM is to close the gaps 



in the rings of the SRRs. This destroys their magnetic 
resonance without substantially affecting their electric re- 
sponse. Closing the gaps, i.e. eliminating the capacitors 
from the rings of the SRRs, the induced circulating cur- 
rents inside the SRRs are still allowed to flow but cannot 
oscillate independently of the external electromagnetic 
(EM) field anymore. The changes in the electric response 
from the closing of the SRRs are expected to be weak (in- 
deed luq moves at most by 3% as can be seen in Fig. 1(c)), 
since only a very small amount of metal is added (in the 
gaps of the SRR) . We have checked all these ideas by cal- 
culating the transmission with closed-SRRs. The results 
(see the dashed lines in Fig. 1) are almost the same as the 
ones for the open SRRs. The only important difference 
is the disappearance of the first dip in Fig. 1(a), and the 
first peak in Figs 1(b) and 1(c), as well as of the peak at 
about 0.09. The features that disappeared are due to the 
magnetic response of the SRRs. This was also confirmed 
with our retrieval procedure, which clearly showed that 
indeed the dip in Fig. 1(a) and the peaks in Figs 1(b) and 
1(c) are due to the magnetic response of the SRRs. The 
most important point of this analysis is that the combi- 
nation of the wires with the SRRs can lower, in principle, 
the uip of the isolated wires to values very close to the fre- 
quency of the magnetic response, ui m , of the SRRs. In 
some cases it is possible that ui' p is lower than the ui m , 
and transmission peaks appear at low frequencies which 
are not LH but RH. This has been also confirmed exper- 
imentally [16]. 

Through our detailed transmission studies and the re- 
trieval procedure for obtaining the effective e and /i of the 
different structures, we were able to obtain analytic ex- 
pressions for the effective e and fj, for all the structures. 
In particular, as expected from effective medium argu- 
ments, an array of metallic wires exhibits the frequency- 
dependent plasmonic form 
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For the frequency dependence of the magnetic response 
of the SRR we obtain in the case of a single resonance 
[17] that 
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We have argued above and we have shown by introduc- 
ing cut-wires instead of SRRs that the electric response 
of the SRRs is equivalent to that of cut-wires. Therefore, 
for the frequency dependence of the electric response of 
the SRR we obtain that [17] 
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(3) 



The magnetic resonance of the SRR is due to the oscil- 
lation of circular currents inside the metallic rings, and is 
determined by the inductance of the loop (enclosed area) 



and the capacitance of the gaps (mainly the gap-width) 
in the rings. The second major contribution of the SRR, 
its electric cut-wire resonance, is due to the oscillation 
of linear currents along the sides of the SRR which are 
parallel to the external electric field. It is basically de- 
termined by the size of the SRR in the direction of the 
electric field. In comparison with the response of the 
continuous wire, Eq. (1), for a cut-wire the electric reso- 
nance is shifted from zero to some ljq ^ because of the 
additional depolarization field due to its finite extent. 
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FIG. 2. Each row shows the response of each one of five 
periodic systems to electromagnetic wave. This response is 
shown through the frequency dependence of transmission co- 
efficient (left panels), dielectric function (middle panels) and 
magnetic permeability (right panels). 

The electric response of the LHM is to first approxima- 
tion equal to the sum of the electric response of the wires 
(given by Eq. (1)) and the electric cut- wire response of 
the SRRs (given by Eq. (3)), which may be approximated 
by the closed-SRRs. The electric interaction between 
wires and SRRs has to be taken into account in order 
to be more realistic. Since the wires or the cut-wires do 
not have any magnetic response, the magnetic response 
of the LHM is just the SRR magnetic resonance, given by 
Eq. (2). In Figure 2 we present in a compact schematic 
form the response of the various "components" of a meta- 
material of SRRs and wires to an EM field. This response 
is shown through the frequency dependence of transmis- 
sion coefficient (left panels), dielectric function (middle 
panels) and magnetic permeability (right panels). The 



first row shows the response of a periodic system of infi- 
nite wires. This response is analogous to that of a bulk 
metal, i.e. there is a cut-off frequency u p above which 
e(u>) becomes positive from negative, and thus the system 
becomes " transparent" to the EM radiation. The system 
does not have any magnetic response. Second row shows 
the response of a system of cut-wires (wires much sorter 
than the wavelength of the EM wave), or closed-SRRs. 
The difference with the continuous wires is that here the 
negative e regime has also a lower edge u>o 7^ 0, due to 
the finite nature of the wires [15]. The e(oj) has the form 
shown in the middle graph. Again no magnetic response. 
Third row shows the response of a periodic system of 
(single-ring) SRRs. Their electric response is cut-wire 
like (like in second row) and their magnetic response has 
a resonance at u> m , where the magnetic permeability (//) 
jumps from positive to negative values. The transmis- 
sion (left panel) becomes finite in the regions of positive 
product e/i and goes to zero for negative e/i. The relative 
order of the characteristic frequencies (w m , uj p , ujo) de- 
pends on the parameters of each specific system. Fourth 
row shows the response of a system of infinite wires plus 
cut-wires (closed-SRR). This system contains in fact all 
the electric response of the LHM, without its magnetic 
response. The combined electric response of wires and 
cut-wires leads to a new cut-off frequency, u>' , much lower 
than u p . Last row shows the full response of the LHM; 
its electric response is that of a system of wires plus cut- 
wires (fourth row) and its magnetic response is that of 
a periodic system of SRRs (third row). Simulations (see 
Fig. 1) and the retrieval procedure have verified all as- 
pects shown schematically in Fig. 2 (for the explanation 
of the peak at 0.09 see Ref. [17]). 

We have systematically studied the transmission prop- 
erties of LHMs composed of SRRs and continuous wires, 
and also of their various components separately. A re- 
trieval procedure was used to obtain the effective param- 
eters for each case. We found that the electric response 
of a LHM is the sum of electric responses of the wires 
and the SRRs. This changes the existing picture where 
the electric response of a LHM was attributed only to 
the wires, and provides a valid scheme for interpreting 
experimentally observed transmission peaks. 
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